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We present hydrodynamical models for Corotating Interaction Regions, which were used by 
Lobel (2007) to model the Discrete Absorption Components in HD 64760. We also discuss our 
failure to model the rotational modulations seen in the same star. 



1 Introduction 

Discrete Absorption Components (DACs) seen in 
ultraviolet lines are one of the clearest indicators 
of large-scale structure in the winds of hot stars. 
Cranmer & Owocki (1996) attributed these DACs 
to Corotating Interaction Regions (CIRs) which are 
formed when slow and fast streams collide in the 
wind. In their model, CIRs are produced by intro- 
ducing a spot on the stellar surface. The spot can 
literally be a spot, but can also simulate the effect 
of non-radial pulsations, or a magnetic field. 

In the present work, we attempt a quantitative 
fit to the DACs seen in the Si iv A 1395 line of 
the BO. 5 la star HD 64760. This paper describes 
the hydrodynamical models used for this fit and dis- 
cusses the problems in fitting the rotational modu- 
lations, which are an additional feature seen in the 
P Cygni profiles of this star. In a companion paper 
(Lobel 2007), the radiative transfer code and the fit- 
ting of the observed DACs is discussed. 



2 Hydrodynamics 

To construct the hydrodynamical models, we use the 
Zeus3D code (Stone & Norman 1992). We solve the 
3-dimcnsional equations of hydrodynamics, limited 
to the equatorial plane. This approximation is jus- 
tified by the high wsini value of HD 64760, which 
suggests that it is seen (nearly) equator-on. 

In our model, we largely follow Cranmer & 
Owocki (1996). We include the radiative accelera- 
tion due to the spectral lines as well as the rota- 
tion of the star and the wind (through angular mo- 
mentum conservation). The energy equation is also 
solved (including radiative cooling and a floor tem- 
perature), but this turned out to be unimportant in 
the present calculations. We further introduce a cir- 
cularly symmetric spot on the stellar surface, with 
its centre on the equator. This spot is specified by 
its brightness A sp and its angular diameter <E> sp . 

Contrary to Cranmer & Owocki, we allow the spot 
to rotate with a velocity v sp that can be different 
from the rotational velocity i> rot . The main justifica- 
tion for this is that the recurrence time of the DAC 



is measured to be 10.3 d while the rotation period 
is only 4.1 d (see Lobel 2007). A similar assumption 
was also made by Kaufer et al. (2007) in their model 
for the Ha variations in HD 64760. 




Figure 1: The colour scale presents the density con- 
trast of the CIR with respect to the 
smooth wind. The velocity vectors plot- 
ted show the differences with respect a 
smooth wind. The model parameters are: 



.4 



0.1, $ sp = 50°, v sp = v rot /2.5 



The effect of a spot on the hydrodynamics of the 
wind is easiest to understand if there is no rota- 
tion. In that case, there is a sector above the spot 
with a different wind (i.e. a different mass-loss rate 
and terminal velocity) compared to the smooth wind 
around it. When the star rotates, these different 
winds collide and a spiral-shaped density enhance- 
ment is formed (see Fig. 1). Trailing this density 
wave is a region with a velocity plateau, which is 
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Figure 2: Contour lines of the effect of a CIR on the mass-loss rate as a function of spot brightness and 
angle. The plusses indicate the models we calculated. The best-fit model is shown with a red 
square. The left panel shows the values measured from the hydrodynamical models, the right 
panel the results from the approximation in Eq. 2. 



mainly responsible for the formation of the DAC 
(Lobel 2007). One should also note that the CIR 
is not a stream of material particles, but is a pat- 
tern in the wind. If we follow trace particles, we see 
that they can cross the whole width of the CIR as 
they move out almost radially. 

We constructed a large grid of hydrodynamical 
models, from which the DACs were then calculated 
using the Wind3D code. The best fit to the Si iv 
A 1395 line of HD 64760 was then determined (see 
Lobcl 2007) . The hydrodynamical model of this best 
fit is shown in Fig. 1. 

The total mass-loss rate of the structured wind 
(M s t r uct) is easy to calculate if we assume the star 
to be non-rotating: 

-^struct = 7-M B p t + ( 1 — — J -^smooth, (1) 

where is the solid angle of the spot, M S mooth is 
the mass-loss rate in the smooth wind and M spot 
that above the spot. Using the approximate rela- 
tion that M oc L 1 /", we can rewrite this as: 

w^*t- a + ^p) 1/a + (i-r)- GO 

M smoo th 4tt V 47ry 

A rotating star will have a structure that is quite dif- 
ferent from the non-rotating model, but close to the 
stellar surface the differences will be small. We can 
therefore use the density and velocity measured near 



the surface to determine the mass-loss rate. The re- 
sults for a series of models are presented in Fig. 2. 
Our best-fit model for the HD 64760 DAC corre- 
sponds to an additional mass loss of only 0.6 %. 



3 Rotational modulations 

Encouraged by the good fit we obtained for the DAC, 
we also tried to fit the "rotational modulations" 
present in the dynamic spectrum. These are the 
banana-shaped, nearly-flat absorption components 
seen in Fig. 3. All our models fail however to reach 
the terminal velocity sufficiently quickly. Our failure 
is surprising in view of the good agreement presented 
by kincmatical models such as those by Owocki et 
al. (1998) and Brown ct al. (2004). 

To investigate this further, we switched to us- 
ing the kincmatical models developed by Brown et 
al. We find that these models can indeed explain 
the observations, provided the v ro t/woo ratio is small 
(~ 0.05). In such a case, the rotational modulation 
does reach the terminal velocity sufficiently quickly 
(dotted white line in Fig. 3). For HD 64760 how- 
ever, we have that Vrot/^oo ~ 0.16 and we obtain 
the solid white line in Fig. 3, which does not fit 
the observations. Because we consider the spot ve- 
locity as a free parameter, we can of course intro- 
duce a slower-moving spot (thereby getting a more 
favourable u sp /t>oo ratio). However, the recurrence 
time scale then becomes longer, stretching out the 
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time-axis, and we again obtain a bad fit which is in- 
distinguishable from the high v ro t /«oo solution (solid 
white line). 
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Figure 3: The dynamic spectrum of the Si iv A 1395 
line in HD 64760. One of the ob- 
served rotational modulations is fitted 
with a kinematical model based on Brown 
et al. (2004). The dashed white line 
has frot/^oo = 0.05, the solid one has 
Vrot/i'oo = 0.16 and is indistinguishable 
from the i> S p/t>oo = 0.05 solution. 



4 Conclusions 

Lobel (2007) and the present paper show that the 
DACs of HD 64760 can be well fitted with hydrody- 



namical models for CIRs. The recurrence time of the 
HD 64760 DACs requires that the spots responsible 
for the CIRs rotate considerably more slowly than 
the stellar surface (by a factor of 2.5). The spots 
cannot be fixed on the surface, thereby excluding 
magnetic fields as a possible cause. As suggested by 
Kaufer et al. (2007), a beat-pattern due to non-radial 
pulsations seems to be the only valid explanation for 
the spot. 

The additional mass-loss rate put into the wind 
is very small and no shocks are formed due to the 
(weak) CIR in HD 64760. The CIRs can there- 
fore not be responsible for the clumping discussed 
at this workshop, nor for the X-ray emission of these 
stars. As Owocki (1998) showed, CIRs are sensitive 
to the amount of instability in the wind. Too large 
an amount of clumping could therefore destroy the 
CIRs responsible for the DACs. As the DACs are 
almost ubiquitous in hot-star winds, this imposes 
important constraints on the amount of clumping. 

Finally, our failure to reproduce the rotational 
modulations seen in HD 64760 suggests that we still 
lack a key component in our understanding of these 
stellar winds. 
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